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We study a twisted vortex bundle where quantized vortices form helices circling around the axis of 
the bundle in a "force-free" configuration. Such a state is created by injecting vortices into rotating 
vortex-free superfluid. Using continuum theory we determine the structure and the relaxation of the 
twisted state. This is confirmed by numerical calculations. We also present experimental evidence 
of the twisted vortex state in superfluid 3 He-B. 



PACS numbers: 67.57.Fg, 47.32.-y, 67.40.Vs 

The equilibrium state of a superfluid under rotation 
consists of an array of quantized vortices, which are paral- 
lel to the rotation axis. Similarly, the equilibrium state of 
a type II superconductor in magnetic field consists of an 
array of flux lines parallel to the field. Here we consider 
twisted vortex states where the vortices have helical con- 
figuration circling a common axis. An example is a vortex 
bundle deformed under torsion. One type of twisted vor- 
tex state appears in a superconductingcurrent-carrying 
wire in parallel external magnetic field 0, 0, ■ The cur- 
rent induces a circular magnetic field, which makes the 
field lines helical, and in order to be force free the flux 
lines take the same conformation. Here we concentrate 
on a new type of twisted state, which can occur even 
when the driving field is a constant. Only this second 
type can appear in charge-neutral superfluids, where the 
rotation is not affected by currents. 

In this letter we demonstrate that the twisted vortex 
state appears spontaneously when vortex lines expand 
into vortex-free rotating superfluid. We present analyti- 
cal results for the twisted state using the continuum the- 
ory of vorticity. In particular, we state the equilibrium 
force-free conditions for a uniformly twisted state, and 
find the equations governing the relaxation of a nonuni- 
form twist. We present numerical simulations for both 
the generation and the relaxation of the twist. We dis- 
cuss the stability of the twisted state, which is limited by 
the helical instability 0, 0] . The results are valid in su- 
perfluids and also in superconductors in the limit of large 
penetration depth and no pinning. Finally, we present ex- 
periments that show evidence of the twisted vortex state 
in superfluid 3 He-B. 

Generation: — We study superfluid in a long cylinder 
that rotates around its axis. We assume that initially 
the system is in metastable state, where no vortex lines 
are present. Then a bunch of vortex loops is created at 
some location. They start to expand along the cylin- 
der. A snapshot from our numerical simulation of such 
propagating vortices is shown in Fig.^ Two striking ob- 
servations can be made. First, the vortices form a front, 
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FIG. 1: The formation of twisted vortex state. The vortices 
have their propagating ends bent to the side wall of the ro- 
tating cylinder. As they expand upwards into the vortex-free 
state, the ends of the vortex lines rotate around the cylinder 
axis. The twist is nonuniform because boundary conditions 
allow it to unwind at the bottom solid wall. The figure gives 
a snapshot (at time t = 25Q~ ) of a numerical simulation of 
23 vortices initially generated near the bottom end (t = 0). 
The parameters are 2-kQ.R 2 / k — 86, a = 0.18 and a' = 0.16 
[g, and R/a = 3 x 10 s , which corresponds to T = 0.4 T c in 
3 He-B at 29 bar pressure. For time evolution see Ref. 0. 



where the ends of the lines bend to the side wall. Second, 
the growing vortex bundle behind the front is twisted, be- 
cause the front rotates at a different speed than vertical 
vortex sections. The existence of the front is deduced 
from simulation and experiment. Here we concentrate 
on the twisted state behind the front. 

The equilibrium state of the superfluid consists of an 
array of rectilinear vortex lines at areal density n v = 
20,/ k, where f2 is the angular velocity and k the circula- 
tion quantum || . The superfluid velocity v s at the loca- 
tion of each vortex is precisely equal to the normal fluid 
velocity v n = fi X r, so that the array rotates rigidly to- 
gether with the cylindrical container. In contrast, the su- 
perfluid velocity vanishes in the vortex-free state, v s = 0. 
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All velocities are given in the laboratory frame. 

Next we determine the rotation of the vortex ends on 
the side wall. For simplicity we consider zero tempera- 
ture, where the velocity of a vortex line i>l is equal to the 
local superfluid velocity. The average superfluid velocity 
in the vortex front is the average of v s between the vortex 
state and the vortex-free state. This gives that the aver- 
age angular velocity of the front is half of the container 
velocity, (</>) = ft/2. Thus the vortex front lags behind 
the vortex array which gives rise to the twisted vorticity. 

It is interesting that the angular velocity of the front 
can also be determined from alternative arguments. One 
is based on the Hamiltonian equations <f> = dH/dL and 
L = —DH/d<j), where L is the component of angular mo- 
mentum along the cylinder axis z. Shifting the vortex 
front vertically, the former equation gives (j> = AE/AL, 
where AE is the energy difference and AL the angular 
momentum difference between the vortex and vortex-free 
states. Evaluating these using the continuum model of 
vorticity (v$ « ft x r in the vortex state) gives the same 
result as above. This result is also easy to generalize to 
the case where the vortex number N is smaller than in 
equilibrium, and the result is 



Nk ln(R/R v ) + 1/4 
~2^ R 2 - R 2 J2 ' 



(1) 



where R 2 = Nn/2irfl. A third argument relies on the 
Josephson relation, where the rotating vortex ends cause 
a phase slippage to compensate the chemical potential 
difference between the two states The rotation of 
one vortex in agreement with Eq. Q has been observed 
experimentally by Zieve et al. in a cylinder with a wire 
on the axis and zero applied flow |10| . 

Uniform twist: — We construct a description of the 
twisted vortex state using the continuum model of vor- 
ticity 0, 0, . We start by considering a twisted state 
which has translation and circular symmetry. This limits 
the superfluid velocity to the form 



i^(r)</> + v z (r)z, 



(2) 



in cylindrical coordinates (r, cj>, z). It is straightforward 
to calculate the vorticity oj — V x v s . The motion of a 
vortex (velocity vi/) is determined by the general equa- 
tion 



as x (> n - fig) - a's x [s x (v n - fi s )] . (3) 



This includes the mutual friction between the vortex lines 
and the normal fluid with coefficients a and a' . Here s is 
a unit vector along a vortex line and v s is the superfluid 
velocity at the vortex core. In continuum theory s = 
u> (the unit vector along a;), and v s — v s + i/V x u> 
differs from the average velocity v s by a term caused 
by the vortex curvature In our case it is a small 
correction but is included for completeness. Here v = 



(k/Att) ln(6/a), b is the vortex spacing, a the core radius. 
Note that only the component perpendicular to s of Eq. 
(0 is relevant since i>l parallel to s is of no consequence. 

We require that vortices do not move radially in the 
twisted state ©• This gives the condition 

, n \ ( v <p i dv 4>\ dv * , v { dv z \ 
ifo-Vi) + + i^—j =0. (4) 

Moreover, this condition implies that all frictional forces 
vanish since the twisted vortices rotate uniformly with 
the cylinder, dl = O x r. The only deviation from solid 
body rotation on the average is swirling v s that every- 
where is parallel to the twisted vortices. We conclude 
that Eqs. (0) and l@J represent a family of stable uni- 
formly twisted states. The wave vector Q = uj^/io z r of 
the twist is an arbitrary function of the radial coordinate, 
Q(r). An explicitly solvable case is obtained by choosing 
a constant Q and neglecting v. 



(ft + Qv )r 
1 + Q 2 r 2 ' 



v z (r) = 



v — Qflr 2 
1 + Q 2 r 2 



(5) 



An important property of the twisted states is the flow 
parallel to the axis, v z (r). Assuming that there is no net 
flow gives an integral condition for v z (r). In the case of 
Eq. © this implies v = (fl/Q)[Q 2 R 2 /ln(l + Q 2 R 2 )-l}. 
The deviation of from fir implies that vortices are 
more compressed in the center and diluted at larger r 
compared to equilibrium rectilinear vortices. 

We note that also the Navier-Stokes equations have a 
stationary solution for uniform swirling flow, but only 
under a pressure gradient along z. The twisted state is 
closely related to the inertia wave in rotating classical 
fluids. Various forms of twisted vorticity as solutions of 
the Euler equation have been studied in the literature 

Nonuniform twist: — Next we construct equations gov- 
erning the relaxation of twisted vortices. Now all com- 
ponents (v r , VrjtjVz) of v s are nonzero and functions of r, 
z, and time t. Here we take into account only first order 
deviations from the rotating equilibrium state. The dy- 
namical equations can be formed using again Eq. © for 
t>L, to obtain equations for the radial and azimuthal co- 
ordinates of vortices. These together with the continuity 
equation and u> = V x v s form a closed set of equations. 
The same set of equations has been derivedpreviously 
starting from the dynamical equation for v s R , [w| . The 
essential result is the dispersion relation 0, \m 
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k 2 )<7 
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-ia(l3 2 + 2k 2 m ) 



±iy/a 2 /3 4 - 4(1 - a') 2 k 2 (P 2 + k 2 )^ 



m 



(6) 



The waves giving rise to this dispersion are of the 
form v r = ckJ\((3r) exp(ifcz — iat) and v z = ic[3Jo((3r) 
exp(ikz — iat). Jq and J± are Bessel functions, while 
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77i = l + ^fc 2 /2ftand?72 = 1 + j/(/3 2 + fc 2 )/20. The bound- 
ary condition v r (R) = leads to (3 — 3.83/i?, 7.01/i?, . . .. 
We study the case that the square root in Eq. © is real. 
Then the negative sign corresponds to radial motion of 
vortices (which induces also azimuthal motion). Here we 
are interested in the positive sign, which corresponds to 
twisting the vortex state. In this case the frequency a 
vanishes for vanishing wave vector k. This is in agree- 
ment with our preceding analysis that there is no re- 
laxation of the uniformly twisted state. For a nonuni- 
form twist we have to consider a finite k, but still assume 
k <C R~ l , The dispersion relation JBJ then simplifies 
to (7 = —ik 2 (2fl/(3 2 + is)/d, where d is another mutual 
friction coefficient d = a/[(l — a') 2 + a 2 ]. This limit 
corresponds to the diffusion equation 



dt dz 2 ' 



i (2fi 

D = -A,w 



(7) 



with effective diffusion constant D. Here f(z,t) can be 
either v r or v z . We see that the diffusion gets faster to- 
wards lower temperatures, where the friction coefficients 
approach zero. 

Simulation: — We have tested the previous theory with 
numerical calculations. In the initial state we have placed 
a number of vortices at one end of a rotating cylinder so 
that they bend from the bottom to the side wall. The 
dynamics is determined by_calculating v s in Eq. @ from 
the Biot-Savart integral [15|. We assume v n — ft x r. An 
illustrative case with a small number of vortices is shown 
in Fig. ^ Another case with more vortices is examined 
in Fig. Altogether with the averaged axial and azimuthal 
velocity profiles. 

The essential features in Fig.|2are the vortex free state 
in the upper part of the cylinder, the propagating vortex 
front, and the twisted vortex state that is left behind. 
In the front increases rapidly so that the azimuthal 
counterflow — fir is strongly reduced in absolute value. 
While the vortex front progresses, the vortex ends rotate 
at a lower speed than the cylinder. This generates the 
twisted vortex state. A clear signature of the twist is the 
axial velocity v z . It is downwards in the center and up- 
wards in the periphery (corresponding to a left handed 
twist, Q < 0). At the bottom wall the boundary condi- 
tion prohibits any axial flow. This implies that the twist 
vanishes there. We assume that the vortex ends can slide 
with respect to the bottom wall. Thus the winding gen- 
erated by the front is unwound at the bottom. 

The calculations indicate that the vortex front deviates 
from any equilibrium configuration and probably cannot 
be described by simple analytic theory. On the other 
hand, the relaxing twist seems to obey the diffusion equa- 
tion (J7J). The profile of v z in Fig.|2]can be understood as 
relaxation towards a steady state where v z is linear in z. 

The twist implies superflow parallel to the vortex lines. 
In such a case it is expected that individual vortices can 
become unstable against helical distortion. A calculation 
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FIG. 2: A snapshot of the numerically calculated vortex state 
expanding along z: (a) vortex configuration, (b) axial v z , and 
(c) azimuthal components of the superflow velocity. The 
velocities, plotted as functions of z, are averaged over the 
azimuthal angle and are shown at radii r — nR/6 with integer 
n. Note that vj, — fir changes sign in (c) close to the center 
of the bundle, as predicted by Eq. ||KJ. The simulation was 
started with 203 vortices and the picture was taken after a 
time interval of 6057" 1 . The parameters are the same as in 
Fig.[T]except 2irflR 2 /K = 214, R/a = 1.5 x 10 5 . For clarity, 
r and z coordinates have different scales in (a). 



in Ref. 4] predicts this helical instability to take place for 
rectilinear vortices when the velocity of the parallel flow 
reaches v z = 2\j2flv. The simulations indicate that the 
maximum axial velocity v z (see Fig. |2J remains smaller 
than this limit. It appears that if any tighter twist is 
created in the front, it is immediately relaxed by insta- 
bilities and subsequent vortex reconnections. Note that 
the same helical instability is responsible for flux flow in 
the force- free configurations in superconductors 0, ■ 

Experiment: — We now turn to the evidence for the 
twisted vortex state from NMR measurements on a rotat- 
ing sample of 3 He-B. The experimental details have been 
described in Ref. . What is essential is that the NMR 
absorption is measured at two symmetric locations near 
the ends of the long sample cylinder. A measuring run is 
shown in Fig. [3] The vortices are injected in the middle 
of the cylinder at time t = 0. Then the NMR line shapes 
change in both locations simultaneously from the initial 
vortex-free form (N = spectrum in the inset of Fig. |2J 
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to that of the final equilibrium vortex state (N 
spectrum). During the transition the absorption is first 
shifted from a "counterflow" peak to an overshoot in the 
"Larmor region" and later redistributed more evenly over 
the entire spectrum. By tuning one spectrometer on the 
counterflow peak and the other on the peak in the Lar- 
mor region, the timing of the two peaks is resolved in the 
main panel. We see that after a flight time of 22 s, the 
vortex fronts reach the spectrometers. The spectrometer 
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FIG. 3: NMR absorption signals as a function of time af- 
ter injection of vortices at t — 0. The two NMR absorption 
records in the main panel show the reduction in the coun- 
terflow peak and the overshoot in the Larmor region. The 
former is interpreted as the arrival of the vortex front and the 
rapid increase in v$. The latter arises from the axial flow ve- 
locity v z , caused by twisted vortices, and the subsequent slow 
relaxation towards the equilibrium state. The insert shows 
the NMR line shapes in the initial vortex-free state iV = 
and the final equilibrium vortex state = N eq . The spec- 
tra are measured at constant temperature and have the same 
integrated absorption. 



tuned to the counterflow peak sees a rapid drop in ab- 
sorption. The other spectrometer tuned to the Larmor 
peak records first a rapid increase, followed by a slow ex- 
ponential relaxation with a time constant of 14 s towards 
the final level of the equilibrium vortex state. 

The quantitative interpretation of the two signals re- 
quires detailed analysis of order parameter textures in 
3 Hc-B with curved vortices and will be presented else- 
where. Very simply though, a large counterflow peak 
comes from counterflow v n — v s that is perpendicular to 
the rotation axis z (when the static magnetic field is along 
z). Conversely large absorption near the Larmor fre- 
quency comes from counterflow parallel to z. Because z 
is a symmetry direction, there is always some absorption 
near the Larmor frequency but it is modest under normal 
circumstances, where no axial flow is present (N = N cq 
spectrum in Fig. EJ. 

The transient absorption traces in Fig. [31 can now be 
understood as a measurement of the vortex state in Fig. [2 
at a fixed detector location 2<jet' The arrival of the vortex 
front at Zdct causes an abrupt reduction in the counter- 
flow peak as \v$ — f2r| is reduced. Simultaneously v z is 
increasing which is seen as a rise in the Larmor absorp- 
tion. The subsequent decrease of v z after the passage of 
the front is seen as relaxation of the Larmor absorption 
towards the equilibrium vortex state. 

The time constant for the decay of the Larmor absorp- 
tion is plotted in Fig. 0] as a function of temperature to- 
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FIG. 4: The time constant for the decay of the Larmor ab- 
sorption (data points) compared to the decay of the twisted 
state according to the slowest mode of Eqs. @ or (line). 
The line has k = tt/£ as fitting parameter giving £ = 18 mm. 
For comparison, the half-length of the cylinder is 55 mm and 
the detector coil is placed between 5 and 15 mm distance from 
the end plate [16J. 

gether with the slowest mode from the diffusion equation 
0. We note that the theoretical eigenvalue is in order 
of magnitude agreement with the measurements. It is 
especially noteworthy that relaxation gets faster with de- 
creasing temperature. This may at first seem surprising 
since the relaxation is usually associated with friction, 
which decreases with decreasing temperature. Our sim- 
ulations, which are time consuming at the experimental 
parameter values, yield a time constant which is larger 
but within a factor of 3 from the experimental value. 

In conclusion, in rotating superfluid a front followed by 
a twisted vortex bundle is the preferred configuration of 
vortex expansion, rather than a turbulent tangle, when 
the induced supercurrents remain below the helical insta- 
bility limit. It appears feasible that a similar state can 
be generated in a superconductor in the superclean limit. 
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